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An alternative to the effective field theory approach to treat ghosts in higher derivative theories is 
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degrees of freedom that plague, among other theories, the higher derivative gravity models that have 
been proposed to explain cosmic acceleration. We consider the extension of this idea to treating 
a class of terms with order six derivatives, and find that for a general term the Euclidean path 
integral approach works in the most trivial background, Minkowski. Moreover we see that even in 
de Sitter background, despite some difficulties, it is possible to define a probability distribution for 
tensorial perturbations of the metric. 
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I. INTRODUCTION 



The suggestion that the accelerated expansion of the universe may be explained by an 
infrared modification of gravity has fueled renewed interest in higher derivative theories and 
their associated pathologies [l|-l9|. Corrections to the Einstein-Hilbert action built with 
contractions of powers of the Riemann tensor may contain four or more time derivatives 
acting on the physical field, the metric. A few special cases aside, systems with more than 
two time derivatives can be described via a classically equivalent Lagrangian that quite 
generically contains ghosts - degrees of freedom with the wrong-sign kinetic terms - leading 
to catastrophic instabilities if they appear in the perturbative spectrum. Unless a scheme 
to deal with the ghosts is chosen, such models are therefore unsuitable to describe physical 
phenomena. 

The best known and standard way to make sense of theories with higher derivatives is 
through the effective field theory approach. From this point of view, terms that contain 
higher derivatives appear from an expansion of an unknown UV-complete theory. This 
expansion, by definition, is supposed to provide an accurate description of the full theory 
only at low energies, and the physical degrees of freedom are assumed to be only those that 
appear in the ground state of the theory [lO|. Classically, the presence of extra solutions to 
the field equations that correspond to the existence of ghosts is considered as an artifact of 
the effective theory, due to the truncation of an infinite series. If then, for instance, it is 
possible to push the masses of these degrees of freedom beyond the cutoff - the energy scale 
below which the effective field theory is trusted - the ghosts can be ignored. 

Despite the ubiquity of the effective field theory idea, alternative procedures have been 
proposed to deal with higher derivative terms in the action 12|. In this paper we will 
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who demonstrated that the 



focus on the prescription introduced by Hawking and Hertog 
Euclidean path integral formulation of the quantum theory allows one to define a probability 
distribution for a scalar field that appears in the Lagrangian with four time derivatives. 

The theoretical differences between these ways of treating theories with ghosts are inter- 
esting in their own rights. However, it is important to note that if the higher order terms are 
considered as corrections to the second order action for the field, the results calculated in 
the Euclidean path integral approach could lead, in principle, to a different physical result. 
In fact, corrections to the probability for the fields may have a different dependence on the 
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"coupling constant" than the equivalent corrections calculated via the effective field theory. 
Here by "coupling constant" we mean the parameter that controls the strength of the higher 
order term in the action. For example, in higher order theories of gravity the "coupling con- 
stant" contains appropriate inverse powers of the cutoff scale, and a different dependence 
on the behavior of such corrections may shift the energy at which they become important. 
This, at least in principle, holds out the hope of an observational test of these competing 
ideas. However, at fourth order, the analysis in jlS] of a term proportional to the Weyl 
tensor squared has demonstrated that in a de Sitter background there is no discrepancy 
between the Euclidean path integral procedure and effective field theory one. 

In this paper, we explore further whether the Euclidean path integral approach can 
be extended to apply generally, and to explore whether observational differences from the 
effective field theory approach can be realized in practice. Since the study of a general 
higher derivative correction to the Einstein-Hilbert action is prohibitively complicated, we 
therefore focus here on a nontrivial correction beyond 4th order, namely the 6th order term 
Vq.-R^i/V^-R'^''. As we will see, this is sufficient to draw interesting conclusions. 

We demonstrate how, in principle, to apply the Euclidean path integral prescription to 
sixth order terms. However, the question of whether it can be applied to a specific system 
such as General Relativity (GR) plus fourth and sixth order corrections is highly-dependent 
on the choice of background. In particular, a Minkowski background always admits choices 
for the "coupling constants" that yield a well defined Euclidean theory, while a de Sitter 
background, due to its explicit time dependence, introduces some complications, since the 
simple requirements to apply the prescription are not met. Nevertheless, as happened in the 
fourth order case, we shall see that this does not preclude the possibility of finding a viable 
result. 

The paper is organized as follows. Section |IT1 is devoted to a brief review of the Euclidean 
path integral approach, and a discussion of its generalization to sixth order for a certain class 
of quadratic Lagrangians. In section UTTl we derive the perturbed action for tensorial modes 
coming from a sixth order action, about two backgrounds, Minkowski and de Sitter. We 
then solve for the classical solutions and perform the canonical procedure to build the path 
integral in the Lagrangian formulation. Finally, in section |IV] we comment on the results 
and present our conclusions. Throughout the paper we use t and r] to denote cosmological 
and conformal times respectively, and denote the time derivative with respect to them with 
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an overdot d/dt = ( ), with the difference between t and r] being clear from the context. 
After a Wick rotation the time coordinate is described by a real parameter that for both 
cosmic and conformal times we call r, and the derivative with respect to it is represented 
by a prime sign d/dr = ( )'. Conformal time is only used when the de Sitter background is 
taken into consideration. Greek indices run from to 3 and Latin indices run from 1 to 3. 



II. REVIEW OF THE HAWKING-HERTOG FORMALISM 

We begin by reviewing the idea behind the Euclidean path integral procedure and dis- 
cussing the ways in which the fourth order case differs from the usual second order treatment. 
In a second order theory the propagator for a field (p defined by a Lagrangian can be 
found computing a path integral between the initial and final configurations 

{{^f;tf)m-u)) = r^dim^wi^sm , (1) 
J (pi 

where the action for the field is given by 

S[<P]= r dt'L{^,<P,t') . (2) 
Jti 

Here 0x represents the state of the field at time tx- A system described by a quadratic 
Lagrangian with a higher number of time derivatives can be transformed into a second order 
system via nonlinear transformations^; for instance a fourth order system with Lagrangian 



can be recast as 



where ipi and ip2 are defined via 



y 777.2 ~ '"^1 V ^^2 ~ "^1 

In the canonical treatment of higher order systems, this transformed Lagrangian (jl]) is the 
starting point and the system is viewed as a multi-field one, where at least one of the newly 



^ We ignore spatial dependence for the moment or, equivalently, we think of the field (/) as a particular 
Fourier mode. 
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defined second order fields is a ghost. In the case at hand it is easy to note that ijj2 has the 
wrong sign for the kinetic term, playing the role of the ghost field. 
The propagator is given by a path integral over both fields 

/'{V'2,-0l)/ 

{iiJ2f,i^lf■,tf)\{i:2^i^u■U)) = / d[tP2{t)]d[iPi{t)]exp[iS[iP2,H] ■ (6) 

Note that via the definitions of ipi and 1^2, this functional integration can be interpreted as 
integrating over the original field and its second time derivative. However, as pointed out 
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12|, this choice presents a problem. For second order systems the propagator obeys the 



composition law 

G(03,0l) = j rf[02]G(03,02)G(02,0l) , (7) 

where G(03,0i) is the propagator between the two states "1" and "3", and "2" represents 
an intermediate state. When one joins the fields above and below the intermediate time t2, 
the value of the field 0(^2) is fixed, but its time derivative is not, resulting in a jump in (l>(t2), 
which in turn corresponds to a delta function in the value of 0(^2)- Unfortunately, in the 
original fourth order action ([3]) the second time derivative appears quadratically, and hence 
the original composition law of the path integral is lost and infinities arise. This argument 
applies quite generally to higher derivative systems. In fact, the standard way to deal with 
fourth order systems is to use Ostrogradski's theorem to define a Hamiltonian from the 
fourth order Lagrangian and to take and to be the canonical variables over which one 
integrates in the path integral. 

On the other hand, in the alternative procedure proposed in to deal with fourth 
order systems, the fundamental variables are taken to be the field and its first time deriva- 
tive. This choice is motivated by the need to retain the continuity properties of the path 
integral formulation, as described above (see 3] for details). However, this point of view 
introduces a different problem. Initial and final states are then described in terms of and 
0, which behave much like position and momentum for a particle in quantum mechanics. 
The proposed procedure is to rotate the system to Euclidean time, and then to integrate 
out the dependence in the definition of probabilities, thus obtaining well defined quantum 
mechanical observables at the price of a loss of unitarity. This procedure has always been 
possible in the special cases studied in the literature so far. 

Therefore, a summary of the practical procedure is: 
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1. From the fourth order action 5* perform a Wick rotation to obtain the Euchdean action 



2. Derive the Euchdean equations of motion and corresponding solutions, 

3. Use the Euchdean version of the path integral to find the propagator for with 
boundary conditions on (f) and </>', 

4. Define a "wavefunctional" as the propagator from a vacuum state at minus infinity in 
Euclidean time, 

5. Find the modulus squared of the wavefunctional, or probability amplitude, which gives 
the probability that a quantum fluctuation leads to a state with specified and 0', 

6. Finally, and crucially, trace over 0' before returning to real time. Note that if one 
were to rotate back to Lorentzian time before integrating, the probability would be ill 
defined, reflecting the existence of the ghost degree of freedom. 

There is no magic trick behind all this, since taming the ghost by integrating over the 
infinities that it introduces happens at the price of a violation of unitarity^. The Euclidean 
formulation of the path integral together with the requirement that the fields die off at 
Euclidean infinity ensures that the fields remain bounded in real time. This is similar to 
using a final boundary condition to remove runaway solutions from systems that would 
otherwise contain them. 

Let us examine this procedure in the specific case of the 4th order system discussed 
earlier. Using t for Lorentzian and r for Euclidean time, rescaling the field 0, the action is 



where is an arbitrary small parameter, the "coupling constant" mentioned earlier. 

After a Wick rotation t — )■ zr the action becomes 




(8) 





^ See the original paper [l2[ for a detailed discussion. 
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so that iS — )■ —5''^. When is positive definite the path integral converges giving a 
well defined Euclidean quantum theory. The resulting equations of motion take the form 
D^cj) = 0, where 



and admit solutions 

(j){T) = Ai sinh(Air) + A2 cosh(Air) + A3 sinh(A2r) + A^ cosh(A2r) . (11) 

The path integral for the propagator from state {(pi, (p'^) at Euclidean time — T to the state 
(02, 02 ) at Euclidean time is then 

((0o,0o;O)|(0t,0t;-2^)) = /^*"'^"^rf[0(r)]exp [-S^M] 

■(0,0) 



e-^"['^-l / ' rf[¥.(r)]exp[-5^M] , (12) 



'(0,0) 

where we have used the decomposition (p = (p^i + ^p, with (p^ the classical solution of the 
Euclidean equations of motion for the appropriate boundary conditions. 

The wavefunctional for a state described by {(pQ, 0q) at time r = is then defined via (fT2|) 

as 

v^o[0o,0o] = lim ((0o,0o;O)|(0t,0t;-2^)) , (13) 

which yields 

^o[0o, <Po] = iVexp + - C^l] ■ (14) 

The values of the coefficients A, B, and C, can be found in the appendix, and is a 
normalization factor found by calculating the path integral over the field ip. There has been 
some debate in the literature about how to actually calculate this normalization function, 



and we refer the interested reader to the very clear article by Zerbini and Di Criscienzo 15| . 
and references therein, for a complete discussion. 
The next step is to define a probability 

P[<Po, <P',] ^ ^0% = N' exp [-2A (0^2 + _ (15) 

As we have already mentioned, this would not provide a well defined probability if rotated 
back to Lorentzian time, since A > and the rotation would introduce (—i)'^ = — 1 in front 
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of the 0Q term. Therefore one rotates back to real time only after integrating over 0', to 
yield as a — )■ 



P[M = W — (1 + ma + . . . ) exp [-m (1 + ma + . . . ) 0o] ^ (16) 



after normalizing the probability density. 

How might this procedure be extended to an arbitrary higher order system with a 
quadratic Lagrangian? Since ultimately we wish to consider higher order terms as cor- 
rections to the propagation of the degrees of freedom of a second order Lagrangian, we seek 
a way to generalize this procedure so that an integration over all the extra degrees of freedom 
is performed in order to obtain the final results. Although some of the original motivations 
presented in jl2| for taking fourth order terms seriously are lost in this approach, this point 
of view is nonetheless consistent with the proposed procedure since it corresponds to trac- 
ing over the unobserved degrees of freedom. Guided by the need for a composition law for 
the path integral, we are led to consider the metric perturbation jij and its derivatives 7^'^, 
and 7j" (rather than 7, 7" and j'^^) as the dynamical degrees of freedom in a sixth order 
Lagrangian for Gravity. The rest of the procedure developed in 12| is then unmodified, and 
in principle the only difficulties that appear should be those associated with the explicit 
calculation of the normalization function for the wavefunctional. 



III. SIXTH ORDER CORRECTIONS 



Since fourth order corrections have already been analyzed in jl3| , we focus here on calcu- 
lating the corrections to the tensor part of the two point function coming from a sixth order 
term. 



A. Expanding the action 

Our goal is to take a convenient contraction of Riemann tensors and their derivatives, 
and to expand it to quadratic order in perturbations about a conformally flat background. 
We will then study the action for the perturbations around two important backgrounds - 
Minkowski space and de Sitter space. 

We focus on one of the simplest covariant terms that contains six time derivatives and is 
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quadratic in metric perturbations, 



(17) 



The total action we start from therefore consists of the Einstein- Hilbert term, a cosmological 
constant, two distinct fourth order contributions and the term above 



-I 



-9 



Ml, 



(18) 



where A = for a Minkowski background, and is nonzero for a de Sitter one. While this 
action is quite general, we shall henceforth ignore the R"^ term; its presence does not affect 
the result as we have explicitly checked, and as one would expect since it merely corresponds 
to an additional massive scalar degree of freedom. This can be seen by changing frame via 
a conformal transformation of the metric. 

Writing the fiat Fricdmann, Robertson- Walker (FRW) metric in terms of conformal time 
77, the perturbed metric is then 

dfj.i' ~ 9 fill "l~ 

- e^' {ri,. + 5l5i^ij) , (19) 

where the scale factor e^^'^'^ is equal to one for Minkowski space and equal to {—Hri)~^ in 
de Sitter space. Since the perturbation •jij is traceless and divergenceless ■ju = dijij = 0, 
the first non zero term in the perturbed action is^ 



drjd' X 



e-'^ [ - {l"'f + + 4p) + - 3j% + 2j%p + 2p^) 



+ 7'( - 4pp''' - 20 p' + Up'p - 48p^) + -f%i + 6^%ip' - p) 
+ ^^-Ap^"" - 2Ap' - %f + 12tp - App'^') 
+ 7^(-8p^-^ + 336p^p + 64p^^p + 56p^ - 48/3 V^^^ - 688p2p2 _ 394^3^/// + 36 (p^^^)^ 

+56pp^ + 104ppp^^^ + 72p6)] - o? iyf - 273 + 73J 
+ [6f (p^ + p) - 675(p^ + p) + f (-12/^ + 72p^p + 12p/^^ - 6p^ - 18p^)] 



+ ^e^1f -75-f(p^ + 2p)] 



(20) 



Note that the background equations have not been used in this derivation. We now specialize 
to the two cases of interest. 



We discuss this expansion in the appendix. 
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B. Minkowski background 



Performing a Wick rotation to imaginary time, and focusing on a Minkowski background, 
for which e^'-^-* = 1, the full sixth order action fEUj) reduces to 



S 



E 
M 



drd X 



(3' 



M2 



(21) 



where, for simplicity, we have omitted the indices on, and the argument of the perturbation 
'jijiri). It is convenient to treat the problem in momentum-space by performing a Fourier 
transform on 7 



(2^ 



E4(^)7|(^)e^ 



ik-x 



(22) 



s=± 



where the polarization tensor satisfies en = k^eij = 0, e*j{k) = eij{—k), and elj{k)el*{k) = 
26'''. 

In order to avoid confusion through notation, we will drop all the unnecessary indices. 
The action for the fc-mode then becomes 



dr 



Ml, 



(-|7"f + 3A;2|7'f - 3A;^|7? + 6A;^|7|') -2«' [Wf - 2k^\i\^ + \^\^) 

+Ml,{\i\^-k^\^\^)] , (23) 



where we have used the notation |7(")p = -^'j-^'j* . Varying this action with respect to 
7* yields the Euclidean equations of motion 



(24) 



with 



- dr^ V + /32 J dr^ 



+ 



d^ 



/32 I dT^ 



- r + 2k 



/32 /32 

Solutions to these equations can easily be written in terms of exponentials as 



(25) 



7cf (^) = ciie^'^ + ci2e-^^^ + Caie^^^ + 0226-^'^ + c^ie^'^ + 0326 



A2T 



3— A2T 



-A3T 



(26) 
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with Ai, A2, and A3 given by 



Ai = A; 



A 



2,3 



(27) 



/32 ^ V /34 

Following the procedure highlighted in the previous section we now define a wavefunc- 
tional that describes the probability amplitude of being in a state characterized by 70, % 
and 7q 



A^exp 



(^oo7o7o + ^oi7o7o + ^027o7o + ^io7o*7o 



+An7o*7o + ^i27o*7o + ^2070*70 + ^2170*70 + ^2270*70'; 



(28) 



The coefficients Aji are functions of the three A,, and we present their explicit forms in the 
appendix. It is, in fact, possible to calculate the normalization factor using Forman's 
theorem 16|. However, since this does not change our result, for simplicity we shall ignore 
the contributions coming from in what follows, until a normalization for the probability 
is needed. 

A probability distribution for 70 can then be defined integrating over 7q and 70 and by 
rotating back to Lorentzian time 



(29) 



where the arrow implies rotating clockwise in the complex plane to Lorentzian time. The 
normalized probability expanded for Mpi ^ 1 then gives 



^[70 



M 



pi 



+ 



exp 



k {2a^ + /3) 



pi 



+ ••• I70P 



(30) 



Interestingly, we have encountered no difficulties in extending the Euclidean path integral 
prescription to our sixth order term in a Minkowski background. This straightforward 
extension suggests that it may be possible to extend the procedure to any system with 2n 
derivatives. 
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C. The de Sitter background 



We now repeat the above calculation in a de Sitter background. As we shall see, the ex- 
plicit time- dependence of the background introduces crucial differences in this case. Setting 
A > and the scale factor to be e'' = {—Hrj)^^, the action in Euclidean time and Fourier 
space becomes 



nE 
•JdS k 



dr 



T 



1112 I //2 

7 +7 



Ml, 



(31) 



Note that if we started without the sixth order term (i.e. set /3 = 0) we would have the 
action presented in isl, which is not positive definite. Nevertheless, the authors of fisj ] 
showed that this does not prevent one from following the Euclidean path integral procedure 
and obtaining a well defined result. We will therefore adopt the same point of view here 
and, although we realize that we are dealing with a non positive definite Euclidean action, 
proceed as planned to see if a meaningful result can be obtained. 

It can also be noted that in principle we could obtain a positive definite action if we 
started from a different form for equation f|T8|) . There, in fact, the signs of and have 
been chosen arbitrarily. If we were to change the signs though, the results presented in 
section IIIIBI would not stand. We choose to keep the sign conventions so that the validity 
of the method is preserved in a Minkowski background. 

Defining, for simplicity, z = —kr, the Euclidean equations of motion become 



Dfjiz) = 



(32) 



with 



where 



d^ 6 d^ 
dz^ z dz^ 



+ 



C2 



-3 + 



dz^ 



z dz-^ 



dz \ z"^ ' z^ 



pi 



/3H 



8 + 24 



AM 



pi 



/3H 



2 y^H^ J 



dz'' 



(33) 

(34) 
(35) 
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Solutions to these equations can be found by factorizing the sixth order differential operator^ 
Dq^, and can be written in terms of exponentials and Bessel functions as 

lci{z) —-^1 [sinh [z) — 2; cosh {z)] + A2 [2;sinh [z) — cosh [z)] 

+>l3-z^ Jai {-iz) + AiZ^Yx^ i-iz) + Aqz^x^ {-iz) + Aqz^Yx^ {-iz) , (36) 

where J and Y are respectively Bessel functions of first and second kind. Recalling that 
z takes values in (0, +00) with +00 being the past infinity boundary, in order to find the 
wavefunctional we need to apply a set of boundary conditions analogous to the one described 
earlier, namely 



7(2;) 7(2;) 70 

^'{z) ^ > z ^ +00 and ^(z) % } z ^ +zo 
i\z) ^ i'{z) ^ 7^' 

The relevant classical solution of the equations of motion is therefore 



(37) 



Bi (1 + z) e-' + ^2^54^ i-tz) + B^zm^' i-iz) , 



(2) 



(38) 



where H^^^ represents the Hankel function of the second kind, and the coefficients Bi contain 
the dependence on zq and on the boundary conditions 70, and 70'. 

To calculate the wavefunction it is sufficient to rewrite the action as 



/TO 
■00 



(39) 



so that on the classical path only the first set of terms survives, with the contribution from 
the integral term being zero. Since we are ultimately interested in integrating over and 
7q it is convenient to collect terms and write the wavefunctional schematically as 



dS 



N exp 



- (^o7o7o + Ai7o7o + A27o7o + ^io7o*7o + ^ii7o*7o 



+ >li27"7o + >l2o7o*7o + >l2i7o*7; + >l227o*7o ) 



■ (40) 



The analytic dependence of the coefficients Aij and D on the parameters a, ^, and Mpi/H 
appearing in the action is somewhat complicated and not very instructive, and so we do not 
display this here. 



For details see the appendix. 
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To make progress analytically we now introduce an approximation scheme, taking a, A 
(if the term is considered) and /3 to be of order unity, with Mpi/H ^ 1 playing the 
role of the small parameter in a series expansion. Beside the reasonable choices for the 
parameters in the action, an extra assumption is needed to simplify the calculation. We 
assume that < 2a^, allowing us to approximate the frequencies Ai and A2 and the Hankel 
functions. With these approximations the associated probability takes a form similar to 
that of equation fHOj) . with the same kinds of terms and different coefficients. In particular, 
focusing on the coefficient of 7o*7o5 which we require to have a negative real part in order 
to proceed with the integration, we find 



P[7o,7o,7o]=AriV*exp 



Pr2 - 1 



2a^ 

/ II* II I 

7o 7o + 



(41) 



P is not yet the probability we are looking for, since integration over 7q and 7q is still needed. 
The bars are a reminder of this fact, counting the maximum number of derivatives acting on 
7o. From equation fHTl) we note that gaussian integration over the real and imaginary parts 
of 7o can be performed only if {kr)"^ > 1. Recalling that k'^f]'^ = k'^/{aH)'^, with a being 
the scale factor, considering fc^r^ > 1 means that the treatment can be considered valid for 
subhorizon modes. 

With the above assumptions both the integrations over 7q and 7q can be performed, and 
after rotating back to Lorentzian time the full final result is reported in the appendix. Before 
we can say we have found a probability for 70, one last check is necessary: the coefficient of 
|7op, in Lorentzian time, has to be negative in order to have a well defined (normalizable) 
probability. We check this by expanding the argument of the exponential as a series in 
Mpi/H, keeping only the leading contribution 



k^ (1 + 2. 



Mi 



pi 



+ 



M 



pi 



7o7o 



(42) 



2(1 + A;V) ' " \ H 

where the symbol L is a reminder that we have rotated back to Lorentzian time. We can 
see that the probability can be integrated over all values of |7o| giving a sensible extension 
of the method in 12| to the sixth order case. This may be compared with the equivalent 



form for the probability in GR, 

PgrIio] 



lA^pexp 



k'M'pi 



2m (1 + Fr/2) 



l7o| 



(43) 
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Finally, from the probability distribution we obtain the two point function for the tensorial 
perturbations 70 in the sixth order case 



The Euclidean path integral prescription is a method to integrate out the infinities appear- 
ing in higher derivative theories with ghosts and extract meaningful probability distributions 
for the non-ghost degrees of freedom. In this paper we have reviewed the original fourth 
order version of the method and have shown how to extend this to a sixth order system 
in a Minkowski background and in a time dependent one - de Sitter. The two cases are 
treated separately since we have shown that a time dependent background, even if highly 
symmetric, introduces some difficulties. The Euclidean action is in fact not positive definite, 
raising doubts about the validity of the underlying quantum theory. Fortunately, as in the 
fourth order case, this does not prevent us from finding a sensible result. 

With higher order gravity in mind, in this paper we have examined an action containing 
GR, a sixth order term and two fourth order ones, with relative strengths set by the Planck 
mass and their relative mass dimension. We have found that the Euclidean path integral 
prescription can be applied to find corrections to the probability distribution of the tensorial 
perturbations about both Minkowski and de Sitter backgrounds. The corrections we have 
found are at least of order one in the de Sitter case, depending on the values of the parameters 
appearing in the action. Therefore the results pose stringent constraints on either the validity 
of the approach, or the presence of the covariant sixth order term considered. 

It is important to be clear about the assumptions made throughout this paper. The first 
one has already been mentioned, and concerns the validity of the quantum theory when the 
Euclidean action is not positive definite. However, note that we could have performed the 
whole calculation in Lorentzian signature, and the present procedure is merely an ad hoc 
prescription for rotating to Euclidean signature only when needed to integrate over ghosts. A 
second problem arises due to the fact that we have chosen 7 as one of our dynamical variables. 
This is somewhat in contrast with the original idea of preserving the continuity properties 
of the path integral. We leave to future studies the analysis of the effect of this particular 




(44) 



IV. CONCLUSIONS 
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choice of dynamical variables. Third, we have considered the simplest possible scheme for 
taking the limit in which the higher order terms become less important in the action; with 
this choice the behaviors of the fourth and sixth order terms are locked together. A general 
approximation scheme in which the two terms may go to zero independently and introduce 
different corrections requires further study. Finally, note that we have only considered one 
specific sixth order term in the covariant action for gravity. Although a full calculation is 
needed, we do not expect the other sixth order terms to conspire and drastically change the 
results found here. 
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Appendix 



A. The Fourth Order Scalar System 



In the fourth order scalar case described by the Euchdean action 



the wavefunctional is defined as 



*o[0o, 0o] = lim ((00, 0o; O)|(0T, 0t; -T)) 

T->-oo 



= A^exp 



Vl — 4m^Q;2 ,2 2m^Q; — tti , my/1 — Aw?a^ 



2(A. - AO '^'^ + ~a{\,-\{)-^'^'' ~ 2a(X,-X,) 
where Ai and A2 are found by solving the equations of motion 



(45) 



, (46) 



Ai = ^ ^{1 - VT- 4m^a^) , A2 = )J ^{1 + VT- Am^a^) . (47) 
The normalized probability, after integrating over 0q and rotating to Lorentzian time, be- 



comes 



P\ 



m Vl — 4m2a2 
ttq; (A2 - Ai) 



exp 



m vT^^4m%^ 2 
a (A2-AO '^^ 



(48) 
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B. Expansion of the 6th order action 



Starting from the action in equation (flSl 



S 



(49) 



the quadratic action for the fluctuations is found by varying the above twice. Since the 
first variation of the volume element is proportional to the trace of 7, and thus zero, the 
remaining terms can be written schematically as follows 



-9 



M2 



■ (50) 



The first few terms, involving the Ricci Scalar, are standard, and the Weyl squared term 
was already calculated in [3], and can be written as 

Finally, the variation of the remaining sixth order term, with the aid of 

S iy^R,,) = h^Si [7/," + 2%- (p - 2p2) - 7,,. + 27., {p''' + 2pp - 4p3) + 2p7i,,fcfc] 



+ 2^a [{^X + i-p'yik + '2iik {p - 2p') + hik,n - (p' - pp)) 



(52) 



gives 



5\V^R,y'^Rn = if""' { - {I'^'f + 7'(-10p' + 4p) + 375 - 375^. + 275(p + 2/3^) 
+ 7'(4p^^ - 4pp^^^ - 20p2 + 44p2^ - 48p^) + 75,, + 675,(p^ - p) 
+ 75(-4p^^ - 24p^ - 8p^ + 72p^p - 4p/^0 

+ 7^(-8p^^ + 336p^p + 64p^^p + 56p^ - 48/9^^^^ - 688p2p2 _ 304p3 
p''' + m{p'''f + 56pp^ + 104ppp^^^ + 72p6)] . (53) 



Combining all of the above, we obtain the full expansion for the action, equation (120 
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C. The Wavefunctional and Probability in a Minkowski Background 



The explicit form for the coefficients appearing in the definition of the wavefunctional, 
equation (1251) . can be cast in terms of Aj as follows 

^0 = A1A2A3 (-2M2«2 + p2 (^_3^2 + ^2 ^ ^2 ^ ^2 ^ X2X3 + A1A2 + A1A3)) 

Aoi = M^i + 2a'Ml (2P + A2A3 + A1A2 + A1A3) + {?>k^ - XlX^ - A^A^ - A? (A2 + A3) 
-X\ (A2 + A3) ' - Ai (A2 + A3) (-3A;2 + A^ + A2A3 + A^) + A2 (3A;2A3 - A^)) 
-20" Ml (Ai + A2 + A3) + (A? + + + A2A2 + XI + XI (A2 + A3) 
-?>e (Ai + A2 + A3) + Ai (A^ + A2A3 + A^)) 
-/32A1A2A3 (Ai + A2 + A3) 
(Ai + A2) (Ai + A3) (A2 + A3) 
A12 = 2Mp>2 - 13^ (-3P + a2 + A^ + A^ + A1A2 + A1A3 + A2A3) 
^20 = /3^AiA2A3 



A02 
Aw 



A21 = (A2A3 + A1A2 + A1A3 
^2 = (Ai + A2 + A3) . 



(54) 



The traced probability in Lorentzian time then reads 

1 



P[7o] = A^A^*exp 



(-4/32AiA2A3(Ai + A2 + A3) (-2M, 



4M2/32(Ai + A2 + A3) 
(-3A;2 + XI + XI + A2A3 + A^ + Ai(A2 + A3))) + {~2M^a\X, + A2 + A3) 

(A? + A^ + A^ + A2A3 + A2A2 + A2A2 + A?A3 + Ai(A2 + A3)' 
-3fc2(Ai + A2 + A3)))')|7or (55) 



D. Classical Solutions and Probability in a de Sitter Background 

The equations of motion in a de Sitter background, equation fl32l) . admit solutions in 
terms of Bessel functions. To find the general solution shown in the text, equation (1361) . 
it is convenient to search for a factorization of the full sixth order differential operator Dq 
defined in equation (1551) . Dq can be split into a fourth order operator acting on a second 
order operator via 

D6[z]^ci{z) = D^[z]D2[z]^ci{z) , (56) 
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where 



1 



1 



2 1 



4^4 



-j 2— h 



25 - 4A^ 
25 - 4A? 



di2 +^ 



25 - 4Ai 



16^6/32/74 



[8Mp^, - _ 48A2 _ 4a2A2) 



+/32i7^(l53 - 104A^ + 16A^)] 



D2[z]=z'^-2zi- 



dz"^ 



dz 



4 



(57) 
(58) 



Here, to avoid confusion, we have replaced the coefficient A^ of the B? term in the action 
with A2, while Aj is a parameter in the decomposition. There are then three independent 
choices of the parameter Aj, namely 



2 ' 



A, 



A. 



13/f2/32 + ^2 


(-4^2-2^ 


1 Ami3^+M^i{Aa^ 


-2/32)-24H2m2j/32(q2+4A2) ^ 






ff2/32 


13/72/32 + Ml 


^-4a2 + 2^ 




-2/32)-24//2M2j/32(a2+4A2) ^ 







1/2/32 



(59) 



(60) 



(61) 



with these choices we obtain the six solutions of f l36|) . 

Once the classical solution is given, it is possible to calculate the associated value of the 
Euclidean action, find a wavefunctional as discussed in the text, and after tracing over the 
unobserved 79 and 79, and rotating back to Lorentzian time, eventually find a probability 
for 7o. 
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The full form of the non-normalized probability is then 



P[7o] = A^exp 



H^k^/3^7]^ HaP (2 + k'^T)^ + kY) r2 



-a^ (1 + kW) rf [--^i^ + + ka^ (1 + kW) r, 

H^/3hirl (18 + 9k'^7]^ + Ak^'^) + Hka/3r] (6 + Yk'^rj'^ + k^rj^) Mpira) 



Ha'^P (8 + 20^)^772 + 21kY + Qk^T]^) r, 



2 



-r — ^ + ka^r^ (1 + A; V) ' 

+-^H^p^rlrl {H^p^ (I8 + Qk'^ri'' + kY) + H^kap'^rj (I8 + 29^7]^ + Wk^ + k^rf) ^pi^2 

+Ha^j3 (1 + A;^?]^ + 5A;^?7'^ + 2^^776) MpVa + ka'^V (4 + Sfc^r^^ + k'^7]^) Mlf l) 

--^Ha^rl (SH^kp^r) (l + k^) + H^oi^^ (l + + - 12A;^?7^ - 4A;^r7^) Mpira 
^pi 

-Hka^/Sf] (-1 + ITfc^r;^ + 21A;V + fc^^^) ^pi^2 + «^ (6 + I7k'^v^ + 18A;V + 7fc^r/6) M^r^) 
+_i_i//3r3^2 (^4^4 ^ 45^2^2 ^ ;l9^4^4^ ^ H^kW/3Y (21 + Qifc^r^^ + A;^r?^) My 



^2 



+2Hka^pr] (12 + nk^r]^ + 5kY + k^) + (l + k^ + k^ + k^) M^^rQ 

(6 + 4A;2ry2 _ 3^4^4) H^ka^^r] (3 + 2A;2ry2 - 2^^774 + k^) 



H^a^p^ (-1 - k^ - 2kY + iOk^ + 4kY) A 
^ Ml 

^H^k(y^,fil (18 + 39A--/r + 26A-i//^ + r>}/'if') ,| ^ 2//A-Vi Khj^ (1 + A-^//^) 



>2 



+A;Q;^r7 (l + k'^r)^) 

+ ^rfr2 [H^P^ (24 + 49^27^2 + iQkY) + H^k^afS^ (I6 + 27A;2r/2 + 9A;V) M^ir^ 
+H^kW/3Y (18 + 9kY + k^) MpVa + H^ka^P^r) (39 + STA;^?^^ + 20A;\^ + Ak^r)^) 
+Ha^P (1 + k^ + Ak^ + 4A;6?7'^) MpVs + ka^Tj (l + A;^^!;^)^ MpVa) 

{H'P'rj (1 + A;2r72) (n - (-//^/3^ (3 + 2kY) n - (3 + 2A;277^) 

+ k^rf) MpVt {Hk''(57f + « (1 + A;'^]') Mpira) 
+Hk^a^Pr)^Myr2 {Hk^Prj^ + a {l + kY^^p,r2) + a^M^ {H^k^^^ 
+3Hk^apr)^ (l + ^7)'') Mp,r2 + a' (l + A;'??')' MpVs) ) ) , (62) 



where ri and r2 are given by 



^M = \/l±\/l-^- (63) 
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